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Measurement of the pion form factor for M£ n between 
0.1 and 0.85 GeV 2 with the KLOE detector 
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Abstract The KLOE experiment at the </>-factory DA$NE has measured the pion form factor in the range 
between 0.1 < M 2 ^ < 0.85 GeV 2 using events taken at y/s = 1 GeV with a photon emitted at large polar angles 
in the initial state. This measurement extends the M 2 „- region covered by KLOE ISR measurements of the 
pion form factor down to the two pion production threshold. The value obtained in this measurement of the 
dipion contribution to the muon anomalous magnetic moment of Aa™ = (478.5±2.0 s tat±4.8 a yst ±2.9 t hoc 
further confirms the discrepancy between the Standard Model evaluation for a M and the experimental value 
measured by the (g-2) collaboration at BNL. 
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1 Introduction 



The anomalous magnetic moment of the muon, 
a M , is one of the best known quantities in parti- 
cle physics. Recent theoretical evaluations [l|, Q 
find a discrepancy of 3 - 4 standard deviations 
from the value obtained from the g-2 experiment at 
Brookhaven 0, @- A large part of the uncertainty 
on the theoretical estimates comes from the leading 
order hadronic contribution a t ^ d ' 1 °, which at low en- 
ergies is not calculable by perturbative QCD, but has 
to be evaluated with a dispersion integral using mea- 
sured hadronic cross sections. The use of initial state 
radiation (ISR) has opened a new way to obtain these 
cross sections at particle factories operating at fixed 
energies 0]. The region below 1 GeV, which is ac- 
cessible with the KLOE experiment in Frascati, is 
dominated by the tt + tt~ final state and contributes 
with ~ 70% to a^ ad ' 10 , and ~ 60% to its uncertainty. 
Therefore, improved precision in the 7T7t cross section 
would result in a reduction of the uncertainty on the 
leading order hadronic contribution to a M , and in turn 
improve the Standard Model prediction for a M . 



2 Measurement of a n7T 

The measurement has been performed with the 
KLOE detector at the DA$NE e + e" collider in Fras- 
cati. DA<I>NE is a ^-factory that usually operates at 
y/s ~ M^, and has delivered ca. 2.5 fb _1 of data 
to the KLOE experiment up to the year 2006, from 
which KLOE has reported two measurements of the 
7T7T cross section between 0.35 and 0.95 GeV 2 0, EJ. 
In addition, about 250 pb _1 of data have been col- 
lected at y/s ~ 1 GeV, 20 MeV below the 4> resonance, 
from which the new results were obtained. Running 
below the (j> resonance diminishes the backgrounds 
from the copious <f> decay products, including scalar 
mesons. As DA$NE was designed to operate at a 
fixed energy around M^, the differential cross section 
dcr(e + e~ — > 7r + 7r~ + 7i S R,)/dAf 2 w is measured, and the 
total cross section = er e +. -_,_+„.- is evaluated us- 



ing the formula ll| : 

d (J ir TT -Vr Ci D 



dM 2 



a_(M 2 J H(Ml,s) , 



(1) 



in which s is the squared e + e center of mass energy, 
and H is a radiator function obtained from theory de- 
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scribing the photon emission in the initial state. Final 
State Radiation (FSR) terms are neglected in Eq. [TJ 
but are taken into account properly in the analysis. 
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Fig. 1. Schematic view of the KLOE detector 
with selection regions. 

The KLOE detector (Fig. [T]) consists of a high reso- 
lution drift chamber {<J p /p < 0.4%) Q and an elec- 
tromagnetic calorimeter with excellent time (a t ~ 
54 ps/^E [GeV] ®100 ps) and good energy (a E /E~ 
5.7%/ y/E [GeV]) resolution 

2.1 Event selection 

The previous KLOE analyses 0, @] used selection 
cuts in which photons are emitted within a cone of 
f? 7 < 15° around the beamline (narrow cones in Fig.[T]) 
and the two charged pion tracks have 50° <6 n < 130° 
(wide cones in Fig. |T|) . In this configuration, the pho- 
ton is not explicitly detected, its direction is recon- 
structed from the tracks' momenta by closing kine- 
matics: p 1 ~p miss = — {p n + +p w ~). While these cuts 
guarantee a high statistics for ISR signal events, and 
a reduced contamination from the resonant process 
e + e~ — > 7r + 7T~7r° in which the n° mimics the miss- 
ing momentum of the photon(s) and from the final 
state radiation process e + e~ — * 7r + 7r _ 7 FSR , a highly 
energetic photon emitted at small angle forces the pi- 
ons also to be at small angles (and thus outside the 
selection cuts), resulting in a kinematical suppression 
of events with M 2 < 0.35 GeV 2 . To access the two 

7T7T 

pion threshold, a new analysis is performed requir- 
ing events that are selected to have a photon at large 
polar angles between 50° < 6 1 < 130° (wide cones 



in Fig. [!}, in the same angular region as the pions 
to be included. The drawback using such acceptance 
cuts is a reduction in statistics of about a factor 5, 
as well as an increase of events with final state radi- 
ation and from <j) radiative decays compared to the 
small angle photon acceptance criterion. The uncer- 
tainty on the model dependence of the <fi radiative 
decays to the scalars / O (980) and / (600) together 
with <j) — > pit — > (7T7)7r has a strong impact on the 
measurement 12J. As an obvious way out of this 
dilemma, the present analysis uses the data taken by 
the KLOE experiment in 2006 at a value of ^/s = 1 
GeV, about 5 outside the narrow peak of the <j> res- 
onance (r^ = 4. 26 ±0.04 MeV 0). This reduces the 
effect due to contributions from f j and gn decays of 
the 0-meson to within ±1%. 
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Fig. 2. MC simulation of Af trk vs. M 2 .^. 
7r + 7r~7 and /i + /x~7 events are located around 
and m M respectively, while tt + tt~tt° events 
occupy a region in the upper left of the plot. 
The black lines represent the cuts used in the 
analysis. 



Contaminations from the processes </>— > 7r + 7r _ 7r° and 
e + e~ — * /i + /i~7 are rejected by cuts in the kinematical 
variables trackmas^ and ^| (see Fig.[5]and Fig. 13). A 
particle ID estimator based on calorimeter informa- 
tion and time-of-flight is used to efficiently suppress 
the high rate of radiative Bhabhas. The radiative 
differential cross section is then obtained subtracting 
the residual background events, N bkg , dividing by the 
selection efficiencies, e seJ (M 2 ^), and the integrated lu- 
minosity using the formula 



dcr^ 



■N, 



bkg 



dM 2 , 



AM 2 

7T7T 



(2) 



* The trackmass is defined using conservation of 4- momentum under the hypothesis that the final state consists of two charged 
particles with equal mass Mtrk an d one photon. 

tfl is the three-dimensional angle between the direction of the selected photon and the missing momentum. 
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where the observed events are selected in bins of 
AIJ, = 0.01 GeV 2 . The residual background con- 
tent is found by fitting the M trk spectrum of the 
selected data sample with a superposition of Monte 
Carlo distributions describing the signal and back- 
ground sources. The fit parameters are the fractional 
normalization factors for these Monte Carlo distribu- 
tions, obtained in intervals of M 2 
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Fig. 3. MC simulation of f2-angle vs. M% v - 
7r + 7r~7-events (blue) are distributed at small 
values of Q,, while the 7r + 7r _ 7r° events (red) 
occupy the region below 0.5 GeV 2 at larger 
values of Q. The black line represents the cut 
used in the analysis. 

2.2 Luminosity 

The absolute normalization of the data sample is 
performed by measuring Bhabha events at large an- 
gles (55° < 9 < 125°), with an effective cross section 
of (Tehabha — 430 nb. To obtain the integrated lu- 
minosity, £, the observed number of Bhabha events 
is divided by the effective cross section evaluate d by 
the Monte Carlo generator Babayaga@NL0 [3, flij . 
which includes QED radiative corrections with the 
parton shower algorithm, and which has been inter- 
faced with the KLOE detector simulation. A detailed 
description of the KLOE luminosity measurement can 
be found in [lit . 

2.3 Radiative corrections 

The radiator function H used to get in Eq. [T] 
is obtained from the PHOKHARA Monte Carlo gener- 
ator, which calculates the complete next-to-leading 
order ISR effects 17[. In addition, the cross section 



is corrected for the vacuum polarisation [18] (run- 
ning of a cm ), and the shift between the measured 



value of M 2 W and the squared virtual photon mass 
M 2 » = (M^) 2 for events with photons from final state 
radiation. Again the PHOKHARA generator, which in- 
cludes FSR effects in the pointlike-pions approxima- 
tion, is used to estimate the latter 19( , and a matrix 
relating M 2 ^ to M 2 , by giving the probability for an 
event in a bin of M% n to end up in a bin of M 2 , is 
used to correct the spectrum. 

2.4 Results 

Using Eq. [1] and Eq. [21 one obtains the two-pion 
cross section a^. The squared modulus of the pion 
form factor \F„\ 2 can then be derived using the rela- 
tiorEJ 

l^( s ')| 2 = -^^(s') , (3) 

where s' = (AfJ 3 .) 2 is the squared momentum trans- 
ferred by the virtual photon and /?„ = \J \ — is 
the pion velocity. 

Fig. [J]shows \F n \ 2 as a function of (M° x ) 2 for the 
new KLOE09 measurement and the previous KLOE 
publication, KLOE08. As can be seen from Fig. [5[ 
both measurements are in very good agreement, es- 
pecially above 0.5 GeV 2 . 
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Fig. 4. Pion form factor \F^\ obtained in 
the present (KLOE09) and the previous 
(KLOE08) analysis. KLOE09 data points 
have statistical error attached, the grey band 
gives the statistical and systematic uncer- 
tainty (added in quadrature). Errors on 
KLOE08 points contain the combined statis- 
tical and systematic uncertainty. 



tin addition, the choice of radiative corrections applied to cr w7r and \F n \ 2 may differ betwen the two. We adopt the definition 
used in 

MM 

. |24|| . in which cr w7r is inclusive with respect to final state radiation, and undressed from vacuum polarisation effects; 
while |-Fjr| 2 contains vacuum polarisation effects and final state radiation is removed. 
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Fig. 5. Fractional difference between \F%\ from 
the KLOE08 and the KLOE09 analysis. The 
band in dark grey represents the statistical er- 
ror of the KLOE09 result, the band in lighter 
grey gives the statistical and systematic uncer- 
tainty (added in quadrature) for the KLOE09 
result. Errors on KLOE08 points contain 
the combined statistical and systematic uncer- 
tainty. 

The cross section corrected for the running of a ern 
and inclusive of FSR, a^Mt is used to determine 
Aa™ via a dispersion integral: 



1 



, - <T 1 \ 
4^3 / 7rir(7) 



ds' 



(s')K(s') , (4) 



where the lower and upper bounds are s min = 0.10 
GeV 2 and s max — 0.85 GeV 2 in the present analysis, 
and the kernel function K{s) is described in 21|. We 
obtain a value of 

Aa™ (0.1 -0.85 GeV 2 ) = 

(478.5 ± 2.0 stat ± 4.8 cxp ± 2.9 thco ) • 10" 10 . (5) 



Reconstruction Filter 


negligible 


Background subtraction 


0.5 % 


f + p7r bkg. 


0.4 % 


cut 


0.2 % 


Trackmass cut 


0.5 % 


7r/e-ID 


negligible 


Tracking 


0.3 % 


Trigger 


0.2 % 


Acceptance 


0.4 % 


Unfolding 


negligible 


Software Trigger (L3) 


0.1 % 


Luminosity (0.1 th (B0.3 exp )% 


0.3 % 


Total exp. systematics 


1.0 % 


FSR resummation 


0.3 % 


Vacuum Polarization 


0.1 % 


Rad. function H 


0.5 % 


Total theory systematics 


0.6 % 



Table 1. List of systematic errors on the Aa^ 
evaluation. 



The evaluation of Aa™ in the range between 0.35 and 
0.85 GeV 2 allows to compare the result obtained in 
this new analysis with the previously published result 
by KLOE 0: 



KLOE Analysis 


Aa™ (0.35 - 0.85 GeV 2 ) x lO" 10 


KLOE09 
KLOE08 


376.6±0.9 stat ±2.4 cxp ±2.1 thco 
379.6±0.4 stat ±2.4 cxp ±2.2 thco 



• KLOE09 
iCMD 



0.15 0.2 0.2S 0.3 0.35 



The two values are in good agreement. 
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Fig. 6. Pion form factor \F%\ obtained in the 
present analysis (KLOE09) and results from 
the CMD and SND experiments. KLOE09 
data points have statistical error attached, the 
grey band gives the statistical and systematic 
uncertainty (added in quadrature). Errors on 
CMD2 and SND points contain the combined 
statistical and systematic uncertainty. 
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Fig. 7. Fractional difference between |F 2 | from 
KLOE09 and the CMD and SND experiments. 
The band in dark grey represents the statis- 
tical error of the KLOE09 result, the band in 
lighter grey gives the statistical and system- 
atic uncertainty (added in quadrature) for the 
KLOE09 result. Errors on CMD and SND 
points contain the combined statistical and 
systematic uncertainty. 
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2.5 Comparison with other experiments 

Fig. M and Fig. show the KLOE09 result for 



\F„\ 2 together with results from the CMD-2 00 
and SND 24| experiments in Novosibirsk. While on 
the p-peak and above, the new result confirms the 
KLOE08 result being lower than the Novosibirsk re- 
sults, below the p-peak the three experiments show 
good agreement. 

3 Forward-backward asymmetry 

The interference in the amplitudes for ISR and 
FSR is odd under the exchange tt + «-> n~ . This gives 
rise to a non- vanishing asymmetry of the distributions 
in the polar angle 9 for the pions A common 

way to express this is the forward-backward asymme- 
try _4 FB : 

N 7 ,+ (6>90°)-N 7r +(6<90°) 



-4FB(M 2 J: 



(6) 



-{6 > 90°) + N 7r +(0< 90°) 

This quantity is an ideal tool to test the validity of 
models used in Monte Carlo to describe the pionic fi- 
nal state radiation. In a similar way, radiative decays 
of the </> meson into scalars decaying into 7r + 7r~ con- 
tribute to the asymmetry [19j, [20j . As can be seen in 
Fig. [5] and Fig. [51 this has a large effect on the asym- 
metry going from data taken at ^/s ~ 1 GeV to data 
taken at ^fs = AT^, especially in the energy region 
below the p-meson mass. Outside the ^-resonance, 
the asymmetry is almost completely dominated by 
the pionic final state radiation, while on the peak of 
the resonance, the decays of the (/)-meson to / 7 and 
also pn contribute significantly. A comparison with a 
Monte Ca rlo p rediction using the PHOKHARA event 
generator [25| with a model for 0-decays and param- 
eters from [26| . together with a pointlike- pion de- 
scription for the pionic final state radiation, shows 
a good agreement with the data for both sets of 
data. Qualitatively, the theoretical descriptions used 
to model the different contributions in the simula- 
tion agree well with the data, although at low M 2 ^ 
the data statistics becomes poor and the data asym- 
metry points have large errors. In particular, the 
off-peak data in Fig. [5] shows very good agreement 
above 0.35 GeV 2 with the pointlike-pion description 
for FSR. Further work is in progress to determine the 
impact of different models for FSR on the asymme- 
try, as well as to estimate higher order effects [27$ . 
In future, the larger dataset from 2004-2005, which is 
almost 10 times larger than the data shown in Fig. [51 
may be used to determine with high precision the pa- 
rameters of the <p decay contributions, in combination 



» 7T U 7T U 7 

and the assumption of isospin symmetry. This will 
then in turn allow to perform a precise on-peak mea- 
surement of the pion form factor down to the produc- 
tion threshold, using the full data sample of about 2 
fb _1 accumulated by the KLOE experiment. 
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Fig. 8. Preliminary forward-backward asym- 
metry for KLOE09 data taken at yfs ~ 1 GeV, 
and the corresponding Monte Carlo prediction 
using the PHOKHARA event generator with 
model and parameters from [261 ]. 
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Fig. 9. Preliminary forward-backward asym- 
metry for KLOE08 data taken at y^s = M^, 
and the corresponding Monte Carlo prediction 
using the PHOKHARA eve nt g enerator with 
model and parameters from [26] . 
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4 Conclusions and outlook 

The KLOE experiment has performed a new mea- 
surement of the pion form factor \F W \ 2 in the Af 2 ^ 
range between 0.1 and 0.85 GeV 2 . The result is in 
very good agreement with the previous KLOE re- 
sult, and extends it down to the two-pion threshold. 
Reasonable agreement was found (especially at low 
energies) with the results obtained from the Novosi- 
birsk experiments CMD-2 and SND. The new KLOE 
result further confirms the discrepancy between the 
Standard Model evaluation for a M and the experimen- 
tal value measured by the (g-2) collaboration at BNL. 

A next step at KLOE will be the measurement of 
the pion form factor using a normalization to radia- 
tive muon events in each bin. In this way, many the- 



oretical uncertainties would become negligible, since 
the radiator function, the vacuum polarisation and 
the absolute luminosity would cancel out in the ra- 
tio of 7T7T7 over events to first order. Pions and 
muons are separated and identified using kinemati- 
cal variables (e.g. the aforementioned trackmass vari- 
able) (2^. The analysis is in a very advanced state 
and a systematic precision similar to the one obtained 
in the absolute measurement is expected. 

The forward-backward asymmetry ^4 FB is an im- 
portant tool to test models for pionic final state radi- 
ation and radiative decays of the 4> mesons to scalars. 
A good check on the validity of models and their pa- 
rameters is crucial for precise measurements of the 
pion form factor below 1 GeV using initial state ra- 
diation, especially when running at the energy of 
^fs = M <t> . 
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